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ABSTRACT 

The recent VLBI observation of the Galactic center black hole candidate Sgr A* at 1.3mm shows 
source structure on event-horizon scales. This detection enables a direct comparison of the emis- 
sion region with models of the accretion flow onto the black hole. We present the first results from 
time-dependent radiative transfer of general relativistic MHD simulation data, and compare simulated 
synchrotron images at black hole spin a = 0.9 with the VLBI measurements. After tuning the accre- 
tion rate to match the millimeter flux, we find excellent agreement between predicted and observed 
visibilities, even when viewed face-on (i < 30°). VLBI measurements on 2000-3000km baselines should 
constrain the inclination. The data constrain the accretion rate to be (1.0-2.3) x lO -9 M0yr _1 with 
99% confidence, consistent with but independent of prior estimates derived from spectroscopic and 
polarimetric measurements. Finally, we compute light curves, which show that magnetic turbulence 
can directly produce flaring events with .5 hour rise times, 2-3.5 hour durations and 40-50% flux 
modulation, in agreement with observations of Sgr A* at millimeter wavelengths. 
Subject headings: accretion, accretion disks — black hole physics — radiative transfer — relativity — 
galaxy: center 



1. INTRODUCTION 

Due to its large angular size, the Galactic supermassive 
black hole candidate Sgr A* is a promising laboratory for 
precision black hole astrophysics using very long baseline 
interferometry (VLBI ) at millimeter wavelengths. Pre- 
vious measurements (|Shen et aTl 120051 : iDoeleman et al.l 
120011: iKrichbaum et alj Il998l )at 7 and 3.5mm found 
a small instrinsic source size for Sgr A* (100-200 /xas 
FWHM) , but at those wavelengths interstellar scattering 
is the dominant co ntribution to the observ ed size. Recent 
measurements by IDoeleman et a l. (2008) at 1.3mm are 
the first at short enough wavelengths to avoid contami- 
nation due to interstellar scattering and at long enough 
baselines to achieve event-horizon scale resolution. These 
measurements have been used bv lBroderick et afl (2009) 
to constrain the black hole and accretion flow parame- 
ters in terms of a 1.5D, nonrelativistic, time-s teady radia- 
tively inefficient accretion flow model (RIAF; I Yuan et al.l 
[200311 . 

Sgr A* is also known to exhibit multiwayelength, polar- 
ized flares (e.g., lYusef-Zadeh et al1l2007t iMarrone et al.l 
120081 lEckart et al.ll2008al : iDodds-Eden et al.ll2009f l. This 
activity has been modeled as hotspots orbiting in the 
inner radii (iBroderick fe Loebl 120051), ad ia batically ex- 
pandin g blobs (|Yusef-Zadeh et all 120061 : lEckart et all 
l2008bh and emissio n from a short, mildly relativistic jet 
(|Falcke et al.lll993l ). 

As a complementary approach to the relatively simple 
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models outlined above, we perform relativistic radiative 
transfer on data from a recent 3D, general relativistic 
MHD (GRMHD) simula tion of a black hole accretion 
disk ([Fragile et al.ll2007f h and compare simulated syn- 
chrotron images to the VLBI observations. In addition, 
we produce light curves at 1.3mm as a first step towards 
explaining Sgr A*'s flaring activity in terms of magnetic 
turbulence. 

In contrast to the results from RIAF modeling, we find 
that for the single black hole spin considered (a = 0.9), 
all inclinations provide excellent fits to the VLBI data. 
Visibility profiles from low (fac e-on) inclinations (i < 
30°) resemble "ring" models (e.g- lDoeleman et al.l2008ft . 
whereas those at higher (edge-on) inclinations are more 
Gaussian. Ongoing work with increased array sensitivity 
should be able to distinguish between these two types of 
visibility profiles. 

Detection of a ring-like profile would provide an impor- 
tant constraint not only on the geometry of Sgr A*, but 
also on the necessary complexity of accurate accretion 
flow models. These results also demonstrate the power 
of mm VLBI to distinguish between accretion models, 
as well as the model-dependence of attempts to extract 
black hole parameters. 

2. METHODS 
2.1. Simulation Data 



iFragile et all 1)20071 ) presented global GRMHD simula- 
tions of black hole accretion disks. They used a spherical 
grid with polar Kerr-Schild coordinates, and the main 
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Fig. 1. — Best fit ray traced images and corresponding 2D, scatter broadened and rotated visibilities for i = 10° (left) and i = 70° 
(right). The intensities are scaled linearly to th e maximum of each panel (from blue to red to yellow), and the UV plane locations of the 
VLBI observations from Doclcman ct al. (2008) are overplotted on the visibilities as triangles. The ring of peak brightness at a radius of 
~ 25/^as in the i = 10° image is due to the transi tion between bound an d unbound photons. Despite the complex structure of the i = 70° 
image, the visibility is similar in shape to that in Brodcrick ct al. (2009) Fig. 2. 



untilted simulation had an effective resolution of 128 3 , 
except near the poles which are underresolved. The 
torus was initialized w ith the analytic, time-steady , ax- 
isymmetric solution of iDe Villiers fc Hawlevl (|2003f ) and 
threaded with a weak poloidal magnetic field with min- 
imum Pgas I Pmag = 10 initially. The magnetorotational 
instability (MRI) arose naturally from the initial condi- 
tions, and the disk quickly became fully turbulent. The 
simulation was evolved for ~ 8000M, or ~40 orbits at 
r = 10M in units with G = c = 1. Only data from the 
final 2/3 of the simulation are used in this letter, once 
the disk is fully turbulent and transient effects from the 
initial conditions have died down. 

This simulation evolved an internal energy equation, 
and injected entropy at shocks. Such a formulation 
does not conserve energy, and produces a more slender, 
cooler torus than conservative simulations which capture 
the heat from numeri cal dissipation of magnetic fields 
(jFragile fc Meierll2009f ). 

2.2. Ray Tracing 

We performed relativistic radiative transfer on the sim- 
ulation data via ray tracing. Starting from an observer's 
camera, rays are traced backwards in time assuming they 
are null geodesies (geometric opt ics approximati o n), us - 
ing the public code described in iDexter fc Ag ol (2009). 
In the region where rays intersect the accretion flow, the 
radiative transfer equation is solved alon g the geodesic 
( Brod erickl |2006[ ) in the form given in iFuerst fc Wul 
( 2004) , which then represents a pixel of the image. This 
procedure is repeated for many rays to produce an im- 
age, and at many time steps of the simulation to produce 
time-dependent images (movies). Light curves are com- 
puted by integrating over the individual images. 

To calculate fluid properties at each point on a ray, 
the spacetime coordinates of the geodesic are trans- 
formed from Boyer-Lindquist to the Kerr-Schild coordi- 
nates used in the simulation. Since the accretion flow is 
dynamic, light travel time delays along the geodesic are 
taken into account. Data from the sixteen nearest zone 
centers (eight on the simulation grid over two time steps) 
were interpolated to each point on the geodesic. Between 
levels of resolution near the poles, data from the higher 
resolution layer were averaged to create synthetic lower 
resolution points, which were then interpolated. Little 
emission originates in the underresolved regions of the 
simulation. 



In this work, we consider sy nchrotron emission from 
Sgr A* using the formulae fromfMahadcv an et al.l (|1996l ) 
and neglecting polarization. The simulation provided 
mass density, pressure, velocity and magnetic field in 
code units. These were con verted into cgs uni t s follo wing 
the procedure described in lSchnittman et ail (|2006l ). ex- 
cept that we scaled the torus mass to match the observed 
flux of Sgr A* at 1.3mm. The length- and time-scales 
were set by the black hole mass, taken to be 4 x 10 6 M Q , 
at a distance of 8 kpc. 

The temperature was taken to be thermal, with the 
ion and electron temperatures equal. Many authors have 
used two temperature models of Sgr A*, with weak cou- 
pling between the dynamically important ions and the 
cooler electrons. However, our single temperature model 
fits the observed flux from Sgr A* at reasonable accre- 
tion rates ~ 10 _9 M Q yr _1 and produces a spectrum near 
millimeter wavelengths that is consistent with observa- 
tions. In this model, the gas is optically thin except at 
high inclinations during flares. The gas temperature is 
~ lO 11 ^ in the regions of peak synchrotron emissivity, 
which is consistent with peak el ectron temperatures in 
RIAF models. We also followed iGoldston etaU (|2005t ) 
and produced images with (T e + Ti)/T e — 10, keeping 
the total gas pressure fixed. These models became opti- 
cally thick before reaching the observed 1.3mm flux, and 
were inconsistent with the location of the peak in the 
quiescent spectrum of Sgr A*. 

3. RESULTS 

3.1. Visibility Fitting 

Our fitting procedur e closely mirrors that described in 
iBroderick et al.l (|2009t l. Simulated images were produced 
for a range of inclination angles, time-averaged accretion 
rates, and observer times (sets of simulation time steps). 
The images were then averaged over 10 minute observa- 
tion periods (8 simulati on time steps) t o matc h the in- 
tegration times used in iDoeleman et aLl (|2008f) . padded 
sufficiently to resolve the shortest baselines, and Fourier 
transformed to 2D visibilities. The visibilities were ro- 
tated to a range of different orientations, £, the projec- 
tion of the black hole spin axis on the sky measured E of 
N. Interstellar scattering was then applied by multiplica- 
tion b y an elliptical Gauss ian as described in Fish et ail 



(2009) using the fits from iBower et all ((20061). Sample 
images and scatter-broadened visibilities are shown in 
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Fig. 2. — Grid of best fit reduced effective x 2 values vs. incli- 
nation and sky orientation £. The scale is from x 2 = 1 (white) 
to 2.5 and greater (black). We find excellent fits at low inclina- 
tions, which are roughly independent of sky orientation. At high 
inclinations, ou r resu lts show a similar shape to that found by 
IBroderick et al.1 ([2009). Overplotted is the probability density vs. 



inclination, marginalized over observer time, accretion rate and sky 
orientation. 
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Fig. 3. — Best fit visibility profiles for low and high inclinations, 
plotted along the line in the UV plane correspondi ng to the average 
l ocati on of the long baseline measurements from Doclcman et al. 
(2008). The visibilities naturally divide into two types. At low 
inclinations, the nearly circular shadow leads t o a minimum in 
the vi sibility profile, similar to the ring model in Doclcman et al. 
(2008). At inclinations > 30°, the profiles monotonically decrease 
with baseline length. A detection in place of the current upper 
limit (open circle with arrow) should favor one set of profiles. 
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Fig. 4. — Upper and lower limits to the average accretion rate as a 
function of inclination at 68% (triangles) and 99% (plus signs) con- 
fidence from probability distributions marginalized over observer 
time and sky orientation for the a = 0.9 simulation. The narrow 
allowed regions demonstrate the ability of VLBt observations in 
conjunction with GRMHD models to constrain the accretion rate 
of Sgr A*. Other simulations may produce different ranges of al- 
lowed accretion rates. 

Fig. [T] At low inclinations, the emission region is es- 
sentially a ring around the shadow, which causes a mini- 
mum in the 2D visibilities. For higher inclinations > 30°, 
Doppler beaming causes the emission to be concentrated 
on one side of the black hole. The visibilities are then 
approximately elliptical, with the minor (major) axis cor- 
responding to the major (minor) axis of the emission re- 
gion. 

Visibilities were int erpolated to the detections in 
iDoeleman et al.l (120081). and an effec t ive y 2 was com- 
puted as defined in IBroderick et~aTl (pflol . Grids of 
minimum effective, reduced x 2 are shown in Fig. [2] over 
inclination angle and sky orientation. At inclinations 
i > 30°, o ur results are in at le ast qualitative agree- 
ment with IBroderick et al.l (|2009f ). The orientation is 
well constrained at high inclinations, since the general 
image shape is nearly static and asymmetric, and the 
long baseline VLBI measurements pick out a specific ori- 
entation for the visibility ellipse. The two distinct bands 
of good fits are due to the (approximate) up-down sym- 
metry of the image. However, RfAF models ruled out 
inclinations i < 30°, whereas for our GRMHD models, 
good fits (reduced x 2 ^ 1-2) are possible at all inclina- 
tions. This is especially evident from the curve of p(i) 
vs. i overplotted in Fig. [2] Inclinations i < 20° are 
less probable due to the sinz prior and large variation 
of probability density with observer time. Although the 
best fit x 2 is roughly the same for low inclinations at all 
£, at any given observer time good fits are restricted to 
less than half the range of sky orientations. 

Visibility profiles from the best fits at many inclina- 
tions are shown in Fig. [3l The profiles represent the 2D 
visibility plotted along the line representing the average 
angl e in the UV plane of the long baseline measurements 
from IDoeleman et al.l ((20081). Also plotted are the data 
from that work. At low inclinations, the profile reaches 
a minimum as described above, whereas at i > 30°, the 
profile decreases monotonically with baseline length. Our 
best fit face-on visibilit y profile is almost identi cal to that 
of the ring model from IDoeleman et al.l (|2008| ). 

Our images are strongly peaked at the circular photon 



4 



Dexter, Agol & Fragile 



4 r -1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — r 




n h .... i .... i .... i .... i .... 1 
5 10 15 20 25 

t (hours) 

Fig. 5. — Lig ht curves of the last 2/3 of the untilted GRMHD 
simulation from Fragile ct al. (2007) at inclinations of 10° (top) 
and 70° (bottom). The i = 70° light curve has been shifted down- 
wards by 1.5 Jy for clarity. Both light curves exhibit consecutive 
flares starting at t ~ 2 hours, which are consistent with those ob- 
served from Sgr A* at mm wavelengths. They are more prominent 
at i = 10° due to higher optical depth between the observer and 
the flaring region at high inclinations. The arrows denote the 10 
minute intervals corresponding to the best fit images shown in Fig. 

m 

orbit, as can be seen in Fig. Q] Low inclination (i < 30°) 
images produced from a single time step of the energy 
conserving simulation bv lMcKinnev fe Blandford (2009) 
are similar to that in Fig. [TJ suggesting that our re- 
sults are not due to the nonconservative nature of the 
simulation. It is not surprising that the VLBI observa- 
tions are well fit by face-on disks. In addition to the 
ring model, a delta function intensity profile at the cir- 
cular photon orbit roughly matches the observations, as 
does t he spherical accretion model from lDexter fc Agoll 
(2009) scaled to the mass of Sgr A*. In fact, any model 
sharply peaked at the circular photon orbit should ap- 
proximately match the visibility data. This is because 
the ratio of the zeroth to first maximum in the visibility 
for such a model is roughly the same as that of the ze- 
roth order Bessel function, and agrees with the ratio of 
the visibility amplitude at ~ 3500MA to the total flux. 

If this GRMHD simulation approximates the correct 
model for Sgr A*, we can also constrain the accretion rate 
for the single spin value used. Upper and lower limits for 
68% and 99% confidence from probability distributions 
marginalized over observer time and sky orientation are 
shown in Fig. 3J Doppler beaming leads to larger fluxes 
at high inclinations, and hence a lower accretion rate is 
necessary to match the observed flux. The limits are 
also more stringent at high inclinations, where the flux 
is less variable (see Fig. [5]). Further marginalizing over 
inclination, we constrain the accretion rate to be (1.0- 
2.3)xl0~ 9 Af Q yr- 1 at 99% confidence. This constraint 
on the accretion rate of Sgr A*, while model-dependent, 
only uses the VLBI visibility measurements. 

3.2. Variability 

A major advantage of using GRMHD simulations in 
place of RIAF models is the ability to probe variability 
from the same dynamical model used to produce spec- 
tra and visibilities. Previous work has used hotspots 
(|Broderick fc Loebll2005l iDoeleman et al.ll2009f) to study 



variability, but this introduces extra free parameters, 
complicating attempts to extract black hole parameters. 
The light curves of our Sgr A* models at 1.3mm shown 
in Fig. [5] are the first from simulation data to consider 
synchrotron emission, and the first to include the time 
delay along geodesies. 

The variability on short timescales is slightly more no- 
ticeable for the i = 70° case. Interestingly, though, the 
longer time scale flaring behavior is more pronounced in 
the i = 10° light curve. The power spectra are well 
described by red noise spectra with power law index 
T = 2.4, 1.7 for i = 10°, 70°. These are both steeper than 
the ob served power law index of T = 1 inlMaucrha n et all 
(2005) at 3mm wavelength. 

The twin flares occur simultaneously in both light 
curves, rise over half an hour and last 2-3.5 hours. The 
flux modulation is 50% (40%) at i = 10° (70°), mea- 
sured from the peak of the flare to the average of the 
light curve immediately preceding it. All of these fea- 
ture s are consistent with mm flare observations of Sgr 
A* (lEckart et alJl2008btlMarrone et al.ll2008tlZhao et all 
120031 : iLi et al.ll2009f ). Since the flares are seen at low in- 
clinations, they are not caused by Doppler shifts from 
hotspots. The peak intensity is attenuated at large incli- 
nation due to rising optical depth. The flares are caused 
by a rise in magnetic field strength near the midplane in 
the inner radii (r ~ 2-6M). They are not due to heat 
from magnetic dissipation since this is not possible in 
the simulation used here. Thus heating from magnetic 
reconnection is not necessary to produce the mm flares 
of Sgr A*. 

4. DISCUSSION 

In this letter, simulated images of a single tempera- 
ture synchr otron emission mode l from 3D, GRMHD sim- 
ulations by iFragile et al.l (|2007| ) have been fit to recent 
1.3mm VLBI observations of Sgr A*. For the single 
black hole spin considered, the high inclination results 
are in at least qualitative agreement with those pre- 
se nted for 1.5D, stationa ry, non-relativistic models used 
bv lBroderick et all (|2009f ). However, we find that in con- 
trast to previous work, images at low inclination provide 
excellent fits to the data, which arc relatively insensitive 
to sky orientation, and present a qualitatively different 
visibility profile due to approximate circular symmetry. 

Radial inten sity profiles from face-on images of recent 
RIAF models (|Yuan et al.ll2009l ) appear more extended 
than those from our GRMHD simulations. The mini- 
mum of the visibility profile for these images occurs at 
to o short a baseline to ma tch the VLBI data. The model 
in iBroderick et al.l (|2009h . on the other hand, has in- 
creasing number density up to the event horizon, in stark 
contrast to the results from the simulations, where the 
number density peaks at r ~ 5M, well outside of the 
event horizon at r — 1AM and the marginally stable or- 
bit at r — 2.3M. Lar ge number densi t ies in the inner 
radii of the model from IBroderick et alj (|2009f ) cause the 
"shadow" to be much s maller than expected (jBardeenl 
Il973t iFalcke et alll2000D . weakening the power on small 
scales (large baselines), and preventing a good fit to the 
VLBI data at low inclinations. 

Increased 1.3mm array sensitivity should favor either 
the low or high inclination profiles in Fig. [3] Our pre- 
dicted profiles would also be sensitive to measurements at 
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1.3mm along baselines of length 2000-3000km. A future 
detection favoring a ring-like profile would constitute the 
first measurement of a black hole shadow. 

The visibility analysis constrains the average accre- 
tion rate of the a = 0.9 simulation to be (1.0- 
2.3)xlO _9 A/ Q yr- 1 with 99% confidence. While this 
constraint is model-dependent, it is independent of 
and consistent w ith constraints from spectral fitting 
(lYuan et al.ll2003|) and polarization measurements (lAgoll 
20001: iQuataert fc Gruzinovll2000t iMacauart et alj|2006t 
Marrone et al.l 120061 ; ISharma et al.l 12007). and in good 
agreement with the previous spectral fits o f single tem- 
peratu re GRMHD synchrotron images by iNoble et al.l 

We have also produced the first light curves of syn- 
chrotron emission from GRMHD modeling of Sgr A*. At 
1.3mm, flares arise naturally from magnetic turbulence, 
and are consistent with observations. Unfortunately our 
analysis is poorly suited for studying variability in other 
wavebands such as the radio, NIR or X-ray. In the ra- 
dio, most the emission originates outside of the initial 
torus pressure maximum at r = 25M, while the NIR 
and X-ray luminosity is produced by some combination 
of synchrotron emission from non-thermal electrons and 
Compton scattering of synchrotron photons. Modeling 
of both is left for future work. In addition, polarization 
measurements from flares are an important observational 
clue. Adding this capability to our relativistic radiative 
transfer code is in progress. 

There are many uncertainties in the analysis presented 
in this letter. While GRMHD allows for a self-consistent 
solution of the primitive variables, it is unknown whether 
the MHD approximation is valid for the low densities at 
the Galactic center. In addition, there are a range of 
possible solutions, depending on the energy prescription 
used. Current simulations also neglect the interplay be- 
tween ions and electrons, which should be handled self- 
consistently. 

If Sgr A* is radiatively inefficient, simulations which 
conserve total energy are probably more appropriate. 



We estimate the time-averaged energy dissipation rate 
from numerical losses in the simulation as ~ 10 37 ergs 
s _1 for an accretion rate of 1.6 x lO -9 M0yr . This 
is larger than the bolometric luminosity of Sgr A*, and 
gives an effective radiative efficiency of ~ .1. The ac- 
cretion flow is advection dominated, in that an order of 
magnitude more energy is lost to the black hole than 
to numerical reconnection. Face on synchrotron im- 
ag es from a single timestep of the 3D, conservative code 
in iMcKinnev fc Blandfordl (|2009h were similar to those 
found here, suggesting that our results may not be spe- 
cific to nonconservative simulations. 

Despite the uncertainties, relativistic radiative trans- 
fer is an important tool for connecting state of the art, 
3D, global GRMHD simulations with observations. Our 
analysis shows that ongoing millimeter VLBI observa- 
tions will not only constrain black hole parameters, but 
accretion models as well. Expanding this analysis to the 
full range of spins will help determine the model depen- 
dence of such future constraints. If Sgr A* is at a low 
inclination, mm VLBI may soon provide the first direct 
evidence of a black hole shadow. 

In addition to the visibility fits and light curves con- 
sidered here, the variability of Sgr A* can be probed 
through measurem ents of closure phases along triangles 
of baselines (e.g., iDoeleman et al.l [2009). Identifying 
salient features of these closure phases is a goal of future 
work. If the flares seen here are robust features of mm 
synchrotron light curves of GRMHD simulations, multi- 
wavelength and polarization measurements can be used 
to constrain emission mechanisms and the non-thermal 
electron distribution. 
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